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Abstract 



The quark-meson-coupling model is used to study droplet formation from the 
liquid- gas phase transition in cold asymmetric nuclear matter. The critical 
density and proton fraction for the phase transition are determined in the 
mean field approximation. Droplet properties are calculated in the Thomas- 
Fermi approximation. The electromagnetic field is explicitly included and its 
effects on droplet properties are studied. The results are compared with the 
ones obtained with the NLl parametrization of the non-linear Walecka model. 



PACS: 21.65.+f, 24.10.Jv, 25.75.-q, 12.39.-x KEYWORDS: Nuclear matter, liquid-gas phase 
transition, relativistic models, quark-meson coupling model 



*e-mail: gkrein@ift.unesp.br 

^e-mail: fscldpm@fsc.ufsc.br 

•^e-mail: mnielsen@if.usp.br 

^e-mail: cp@teor.fis.uc.pt 



1 



I. INTRODUCTION 



One of the most important problems in contemporary nuclear physics and astrophysics 
is the determination of the properties of nuclear matter as functions of density, temperature 
and the neutron-proton composition. In fact, neutron-star matter at densities between 
0.03 fm~^ and nuclear matter saturation density consists of neutron-rich nuclei immersed 
in a gas of neutrons. In particular, understanding the transition crust-core in neutron stars 
is essential for explaining a number of properties of these stars To achieve this goal, one 
must study not only the ground and excited states of normal nuclei, but also nuclear states 
of high excitation and far from stability. 

Recently, two of us studied droplet formation from the liquid-gas phase transition in 
cold 10] and hot [Q asymmetric nuclear matter in the context of the non-linear Walecka 
model (NLWM) [Q. In the present paper we employ the quark-meson-coupling (QMC) 
model, originally proposed by Guichon In the QMC model, the mean scalar a and 
vector uj meson fields couple directly to the confined quarks inside the nucleon bags, instead 
to point-like nucleons as in Walecka-type models 0. The saturation of nuclear matter 
in this model is due to the density dependence of the effective aNN coupling, generated 
by the coupling of the mean a field to the quarks. As we shall show in this paper, this 
density dependence of the coupling has important consequences for the phase-transition 
and for the droplet radius and surface energy density. While the QMC model shares many 
similarities with Walecka-type models 0, it however offers new opportunities for studying 
nuclear matter properties. Perhaps one of the most exciting ones is the possibility of using 
the same model to study nuclear phenomena in a large range of densities. With the QMC 
model, we have the opportunity to investigate the density regime where the quarks remain 
confined in the nucleon bags, but the stucture of the nucleons nevertheless changes, as 
became evident from the EMC effect 0. We also expect to use the QMC model at much 
higher densities, where the nucleon bags start loosing their identity and the deconfinement 
transition starts taking place It is therefore important to explore the performance 

of the model in all such density regimes. Another good reason to use the QMC model in 
the confined regime is that it might be of help to fix the low-energy constants of relativistic 



Lagrangians, where nucleon substructure is incorporated through a derivative expansion [10 



Since the original version of the model, there have been several ameliorations and exten- 



sions. These include the treatment of the nucleon center-of-mass motion [|Tl|,|I^, treatment 
of finite nuclei |13- |l5|] , change of the bag constant in medium |jl6|, inclusion of finite tem- 
perature effects ||T^ and the treatment of Fock [|I8| and quark-exchange [^] terms. For a 
list of applications of the model to a variety of nuclear phenomena and earlier studies within 
the QMC model, see Ref. m. 



Within the framework of relativistic models, the liquid-gas phase transition in nuclear 
matter has been investigated at zero and finite temperatures for symmetric and asymmetric 
semi-infinite systems [pTH2^ and for finite systems 0,0]. The present study aims to study 
the liquid-gas phase transition and droplet formation in a vapor system at zero temperature 
in the context of the QMC model and explore differences from the results of Refs. [0,0. 
We include the Coulomb interaction and work in the Thomas-Fermi approximation. We 
determine the conditions for phase coexistence in a binary system by building the binodal 
section of the QMC model at zero temperature. As shown in Refs. [0,11, the optimal nuclear 
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size of a droplet in a neutron gas is determined by a delicate balance between nuclear 
Coulomb and surface energies. The surface energy favors nuclei with a large number of 
nucleons A, while the nuclear Coulomb self-energy favors small nuclei. However, because 
of the density dependence of the effective aNN coupling in the QMC model, the critical 
pressure and proton fraction of the phase transition and the droplet properties turn out to 
be significantly different from the NLl parametrization of the NLWM [Q. We have chosen 
to compare the QMC model results with the ones obtained with the NLl parametrization 
because this parametrization has proven to give an excellent reproduction of the ground- 
state properties of the nuclei in general [|,^]. The values calculated for the QMC model are 
also parameter dependent, but other possible parametrizations produce the same qualitative 
features as the one chosen in this work. We also see that, as expected, the presence of the 
outside neutron gas reduces the surface tension, since as the density of the system increases, 
the inside and outside matter become more and more alike. 

The paper is organized as follows: in Sect. II we briefly summarize the QMC model for 
finite asymmetric systems including the couplings of the p meson and the photon to the 
quarks. In Sect. Ill we discuss the Thomas- Fermi approximation and in Sect. IV we apply 
the model to asymmetric nuclear matter. Finally in Sects. V and VI we give our numerical 
results and conclusions. 

II. THE QMC MODEL FOR FINITE ASYMMETRIC SYSTEMS 

In a nucleus, the motion of the nucleon is relatively slow and the quarks are highly 
relativistic. The internal structure of the nucleon has therefore enough time to adjust to the 
external local meson fields, the motion of the nucleon can be treated as a point-like Dirac 
particle [|T3H T5[| with an effective mass M'^. The effective mass depends on the position 
only through the mean scalar field a. To describe a finite system with different numbers 
of protons and neutrons, it is necessary to consider also the contribution of the p meson. 
Besides, any realistic treatment of nuclear structure also requires that one introduces the 
Coulomb force. Therefore, a possible Lagrangian density for a static system is 

Cqmc = ^[^1 ■ d - M^(a(r)) - gM^ho - ^?p^Kr)7o - ^{l + T.^)A{r)^o]^ 

+ ^[(V6(r))^ + mj6(rf] + l(VA(r))% (1) 

where ip{r), a(r), uj{r), 6(r) and A{r) are respectively the nucleon, meson a and the mean 
values of the time component of uj, p and Coulomb fields in the nucleon rest frame; /2 is 
the third component of the nucleon isospin operator; rricr, and rrip are respectively the 
masses of the a, uj and p fields; g^j and gp are the u — N and p — N coupling constants, which 
are related to the corresponding quark coupling constants as g^ = 3g^, gp = g^. Finally, e 
is the electric charge. 

The effective nucleon mass is calculated in the MIT bag model. Parametrizing the 
sum of the center-of-mass and gluonic fluctuation corrections in the familiar form —z/R, 
where R* is the in-medium bag radius, takes the form 
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= + \r:BW\ (2) 

where VLq = ^x*"^ + {R*m*y is the kinetic energy of the quarks, m*{r) = niq — fi'^cr(r) is the 
effective quark mass, rrig is the bare quark mass, is the quark-cr couphng constant, and 
X* is the in-medium bag eigenvalue determined from the boundary condition 



Jo (a; 



- R*m* , 

Ji{x ) • (3) 



\Qg + R*m*q 

The bag radius is obtained by minimizing with respect to R* 

^ = (4) 

and the bag constant B and the parameter z are fixed to reproduce the free-space nucleon 
mass (Mtv = 939 MeV). In this work we use rrig = 5 MeV and fix the free bag radius at 
Rb = 0.8 fm. The results for B and z are: B^/^ = 170.0 MeV and z = 3.295. Results for 
B and z for other values of the bare quark mass and bag radius can be found in Ref. . 

The variation of the Lagrangian, Eq. (|l]), results in the following equations for a spheri- 
cally symmetric system 

d"^ 2d 

-7^a(r) + --i-t^(r) - mla{r) = -g^C{(T{r))ps{r) , (5) 
dr"^ r dr 

(f 2d 

+ rdr^^^^ ~ " -9u;PB{r) , (6) 

^^^""^ ^ r ^^^"^^ ~ ""'^^"^^ " -yP3(r) , (7) 

4^Mr) + -^A{r) = -epp(r) , (8) 
ar^ r dr 



where 



with 



BM* 

g.C{a) = = 3glS{a) , (9) 



^_ nj2 + miRsinq-i) 

^^^^ " n,{Q,-l)+m*Rs/2 ■ ^'"^ 

Note that the a — N coupling constant g^ is related to the quark-o" coupling g^ through the 
relation 

g^ = ?>glS{Q), (11) 

and therefore 
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Also, in Eqs. @ to (H), ps is the nucleon scalar density in the system with A nucleons 

Ps{r) = {A\ [V^p(r)^p(r) + V^„(r)V^„(r)] \A) , (13) 

ipp and Tpn are the proton and neutron spinors, respectively, and ps and ps are given in terms 
of the proton and neutron vector densities Pp and Pn 



Pb = Pp + Pn 



P3= Pp- Pn, 



with 



Pp = {A\^lir)i;pir)\A), p„ = {A\^lir)Mr)\A) . 



(14) 



(15) 



III. THE THOMAS-FERMI APPROXIMATION 

Since we want to compare results with Refs. we employ the semi-classical Thomas- 
Fermi approximation, instead of solving the Dirac equation for the nucleons. This amounts 
to assuming that the mesonic fields vary slowly enough so that the baryons can be treated as 
moving in locally constant fields at each point of space. The basic quantity in the Thomas- 
Fermi approach is the phase-space distribution function for protons and neutrons 

fi{r,k) = e{kF,{r) - k) , i=p,n (16) 

where kpX''^) is the local Fermi wave number for protons and neutrons. The total energy of 
the system is given by 



i=p,n 



,2 „2 



- iVcof - rnluj^ - {Vby - my - iy Af , (17) 



where /ij(r, k) = y P + M'^ {a{r)) + z/i(r), with 

^p{r) = gMr) + ^b{r) + eA{r) , (18) 

Mr) = gMr) - . (19) 

The thermodynamic potential is defined as 

Q = E- PiN,, (20) 

i=p,n 

where pi is the chemical potential for particles of type i and Np and Nn are, respectively, 
the numbers of protons and neutrons 

Ni= [ £rpi{r), i=p,n , (21) 
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with the proton and neutron vector densities, pp and pn of Eq. (|15[), given exphcitly by 

MO =2/ (^/^(-'^) = ^^^.W- (22) 

Minimizing the thermodynamic potential VL with respect to the local Fermi momentum, 
the following expressions for the proton and neutron chemical potentials are obtained 

Pp = ^I^^Tm*^ + g^u + ^h + eA, (23) 



f^n = ^Jkj.„ + + g^u - , (24) 

which can be used to find kp-{M]^, u, b, A, pf, r). 

The fields that minimize fl satisfy Eqs. (^ to (^), where the scalar density of Eq. (|13|) 
is given explicitly by 

P.«^2E/^^i£^^/.(n*). (25) 



i=p,n 




IV. LIQUID-GAS PHASE TRANSITION IN ASYMMETRIC NUCLEAR 

MATTER 

Phase transitions in binary systems are more complex than in one-component systems 
because two kinds of instabilities can occur. We have therefore two stability conditions. We 
have the condition for mechanical stability, which requires 

(26) 

where P is the pressure and Yp = pp/ ps is the proton fraction. We have also the condition 
for diffusive stability, which implies the inequalities 

These reflect the fact that in a stable system, energy is required to increase the proton 
concentration while the pressure is kept constant. 

In the mean field approximation for infinite nuclear matter, the meson fields are replaced 
by their expectation values and Eqs. (^ to (||) become (omitting the electromagnetic field) 

a = ^C{a)p, , (28) 



ml 
_ 9w 



UJ 



PB , (29) 



P3 , (30) 
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where the sources of the fields are constants and can be related to the nucleon Fermi mo- 



mentum through equations (^51) and (p2l), with the distribution functions given by (16) 
and the fields replaced by their expectation values. 

Under this approximation, the energy density and pressure are given by 



^ = 2E / + Mn\^) (^(k,. - A:) + f4 + ^ + (31) 



i=p,n 



There are six parameters to be determined: g^,, g^, gp, m„, m^ and m^. We take the 
experimental values of = 783 MeV and = 770 MeV, and the other parameters are 
taken from Ref. |T§: gl/An = 3.12, gl/An = 5.31, g^^/An = 6.93 and m^ = 418 MeV. 



The two-phase hquid-gas coexistence is governed by the Gibbs condition 

Hi{pp,Pn,M;f) = fi,{p^,p^,M*^'), i = p,n (33) 
P{p„p^,M*^) = P{p'^,p'^,M*^') . (34) 

We have made use of the geometrical construction [p6l , p7| in order to obtain the chemical 



potentials in the two coexisting phases for each pressure of interest. In Fig. 1 we plot Pp and 
Pn as a function of the proton fraction for P = 0.1 MeV/ fm^ (solid line). For comparison, 
we also show in this figure the results obtained with the linear Walecka model (dotted 
line) and with the NLl non-linear Walecka model |Q, discussed in Refs. (dashed line). 

In order to obtain the binodal section, which gives the pairs of points under the same 
pressure for different proton fractions obtained using the geometrical construction, we have 
used the condition of the diffusive stability, Eq . (pTl) , and solved simultaneously Eqs. (|33D 
and (p^), together with 

^ = Pn, MIj) , (35) 

a' = ^C(a')p.(p;,p;,M-). (36) 

The binodal section is shown in Fig. 2 for the QMC model (solid line). It is divided into 
two branches by the critical point (CP). One branch describes the system in a high-density 
(liquid) phase, while the other branch describes the low-density (gas) phase. The liquid 
(gas) phase appears at the right (left) of CP. 

For finite temperatures, T 7^ 0, the gas phase is also characterized by 1^ 7^ for all values 
of the pressure. However, for T = and for P ^ 0.43 MeV fm~^, the gas phase always has 
1^ = 0. For the sake of completeness, we also show in Fig. 2 the binodal sections for the 
NLl discussed in refs. (dashed line). We see that the QMC has a CP with higher values 
of P and Yp, but that in the liquid phase it is mainly characterized, for a given pressure P, 
by smaller values of Yp than the NLl. 
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V. NUMERICAL RESULTS FOR FINITE SYSTEMS 



The solution for the infinite system gives us the initial and boundary conditions for the 
program which integrates the set of coupled non-linear differential equations to (|^) in 
the Thomas-Fermi approximation. In this work the numerical calculation was carried out 
with the iteration procedure described in Ref. 0, which uses also as input the size of the 
mesh, Rmesh- The size of the mesh determines the size of the droplet and, consequently, 
the chemical potentials and the number of particles in the droplet. Hence, the number of 
particles and the proton fraction within the droplet are output of the program and not fixed 
from the start. The same is true for the liquid and gas proton fractions. In order to obtain 
a certain number of particles inside the droplet, we vary Rmesh for a given initial condition 
until the desired number is obtained as output. Therefore, by fixing the number of the 
particles in the droplet we can study the behavior of some of its properties, such as the 
surface tension, neutron thickness, proton radius, etc., as a function of the central (r=0) 
proton fraction. For more details about this numerical procedure we refer the reader to 
Refs. 

The droplet surface energy per unit area in the small thickness approximation is given 

by 



S 



surf 



dr 






(37) 



In Ref. the results obtained from this expression for the density surface energy were 
parametrized and compared with the liquid drop model results. The conclusion was that, 
albeit approximate, Eq. (^) provides a very good estimate of the surface energy density. 
The proton and neutron radii in spherical geometry, Ri {i = p,n), are defined as 



R' 1 

Pp{r)r^dr = - pp^gRl + pp^i{R 

o 



/3 



Rt 



(38) 



and 



R' 1 

3 



Pn,gRn + Pn,l{R'^ " R 



(39) 



where pi^i and pi^g refer to the liquid and gas density respectively; R' is the value of r such 
that \f{r) — fg\ < 10~^, with / being either a meson field or the baryonic density at r, and 
fg the corresponding gas value. This means that R' is the value of r for which the fields and 
density reach their asymptotic gas values. 

Another important quantity is the thickness of the region at the surface with extra 



neutrons known as neutron skin. The neutron skin thickness is defined as 



e = R. 



R. 



.p. 



(40) 



In order to understand the differences between the QMC and NLl in the Thomas-Fermi 
approximation, we have calculated the proton density of a droplet with total particle number 
A = 40 and proton fraction Np/A = 0.5. The results are shown in Fig. 3 and Table I. The 
solid and dashed lines in Fig. 3 show the results for QMC and NLl, respectively. It is clear 
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from this figure tliat tlie QMC result sliows liigher central density and smaller proton radius. 
This feature is present in all solutions we have obtained. In average the central density is 
~ 0.01 - 0.02 fm-3 higher in the QMC model. 

In Table I we give the central density, the surface energy density, the surface thickness and 
the proton radius of the droplet, as well as some properties of the models used, namely the 
incompressibility, efective nucleon mass and the scalar meson mass. The surface thickness, 
t, is defined as the width of the region where the density drops from 0.9pso to O.lp^o, where 
pBo is the baryonic density at r = 0, after subtracting the background gas density. From 
Table I it is seen that the QMC surface energy density is higher and the surface thickness 
smaller. In fact, as a common trend in all our solutions, the surface thickness is ~ 0.4 — 0.5 
fm lower in the QMC model. This, in part, can be explained by the different values of the 



incompressibility K, of the models Intuitively, the higher the value of K, the sharper 
the surface and more similar to a theta function the density distribution becomes. However, 
as discussed in [^] and also verified in [^, this is not always the case when non-linear 
couplings are used. We should stress that the surface properties depend strongly not only 
on K but also on and the effective mass M* as shown in |3^. In fact a smaller value of 
rUs corresponds to a larger range of the attractive interaction , and thus an increase of the 
surface extension. 

In Fig. 4 we show the surface energy density for a droplet with total particle number 
A = Np + Nn = 20, with (dashed line) and without (solid line) the electromagnetic field. 
For comparison we also show in this figure the result (without the Coulomb field) for the 
NLl (dotted line). As we can see, the differences resulting from different models are much 
larger than the differences due to the inclusion of the electromagnetic field. The density 
surface energy increases with the central proton fraction. In fact for small proton fraction, 
the matter inside the droplet becomes more and more neutron rich and, at the same time, 
the density difference between the matter inside and outside the droplet becomes smaller. 
Were the matter inside and outside the same, of course there would be no surface energy. 



This effect was already observed in earlier studies of neutron star matter [32|. Comparing 



the QMC and NLl results, we see that the density surface energy in the latter is smaller, in 
part due to its lower incompressibility, as discussed before. We also see that the Coulomb 
field, a long range force, does not change much the density surface energy of the droplet, as 
expected. 

In Fig. 5 we show the density surface energy for droplets with A = 20 (dashed line) and 
A = 40 (solid line), both with the inclusion of the electromagnetic field, as a function of 
the central proton fraction. We see that the density surface energy increases with the total 
number of particles in the droplet. This is because the systems we studied, with A = 20, 
are not big enough and, therefore, Ssurf depends on R, the size of the droplet. In general 
Ssurf is identified with the surface energy density associated with an infinite plane surface 
with a finite thickness, with the matter on one side being in the liquid phase and on the 
other in the gas phase. We have checked that the density surface energy does not change 
appreciably in going from systems with A = 40 to A = 60. 

It should be stressed that when the Coulomb field is present the solution Yp{r = 0) = 0.5 
does not exist in our formalism, since the Coulomb field favors an increase of the neutron 
fraction. 

In Fig. 6 we show the proton radius as a function of the central proton fraction for the 
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QMC model with (dashed hne) and without (sohd line) the contribution of the Coulomb 
field, and for the NLl without the Coulomb field (dotted line) for a droplet with A = 20. 
As the equation of state of the NLl model is softer than the QMC, the proton radius is 
larger in the former. It is interesting to notice that for Yp{Q) in the interval (0.4 — 0.45), Rp 
has a minimum. This is due to a competition between the repulsive Coulomb force between 
protons and the attractive force between protons and neutrons. Of course the inclusion of 
the Coulomb field increases the proton radius (due to repulsion). 

The fact that for small values of the liquid proton fraction, Rp is smaller when the 
Coulomb field is included is because for these droplets the total proton fraction in the droplet 
is in fact smaller than the total proton fraction in the droplet without the Coulomb field, as 
can be seen in Fig. 7 where we plot the total proton fraction in the droplet {Ytot — Np/A) 
as a function of the central proton fraction in the system. 

Finally, in Fig. 8 we show the neutron skin thickness as a function of the central proton 
fraction for a droplet with A — for the QMC model with (dashed hne) and without (solid 
line) the contribution of the Coulomb field, and for the NLl model without the Coulomb 
field (dotted line). We see that the behavior is very similar in both models, the neutron skin 
being thicker in the latter because of its softer surface. We also see that the contribution of 
the Coulomb field decreases the neutron skin thickness, since the proton radius Rp increases. 

It is worth pointing out that no results were shown for very small proton fractions 
because, for the desired total number of particles (20 and 40) we did not abtain convergence 
for the coupled differential equations. In fact, for small proton fractions, when the Coulomb 
field is switched off we obtain convergence for sufficiently large droplets. As soon as we 
switch on the Coulomb field the binding energy decreases and no sohitions appear whatever 
the size of the droplet. This occurs in the QMC model, with respect to the NLl, for smaller 
proton fractions. 

VI. CONCLUSIONS 

In this work we have studied, within the framework of the QMC model, the properties 
of droplets which arise from the liquid-gas phase transition under conditions predetermined 
by the binodal section. The droplets are described using the Thomas-Fermi approximation. 
The results are then compared with the NLl. 

Some of the consequences of the electromagnetic field in the droplet formation within the 
QMC model, for a fixed number of particles, is to slightly increase the surface energy density 
and to decrease the central density, increase the proton radius and decrease the neutron skin 
thickness except for very small proton fractions or for almost symmetric nuclear matter. 
Nevertheless, the effect of the electromagnetic field is small in surface properties and the 
differences resulting from the use of different models are by far more important. The most 
important effect of the electromagnetic field is, however, to prohibit the existence of droplets 
with a very small percentage of protons, as well as very large droplets. 

We have seen that the implicit density dependence of the couplings in the QMC model 
has important consequences for the liquid-gas phase transition and the properties of the 
droplets. This, together with the different incompressibihties of the models affect sensibly 
the surface properties of the droplets. In particular, the QMC model predicts a higher 
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density surface energy, a smaller proton radius and a smaller neutron skin thickness for the 
droplets than the corresponding predictions in the NLl. 

We also note that while our numerical results depend on the particular model chosen 
(namely, the QMC or NLl models), some qualitative features such as the increase of the 
surface energy density, the increase of the total proton fraction and the decrease of the 
neutron skin thickness with the central proton fraction in the droplet, apply to both models. 

The inclusion of finite temperature effects in droplet formation within the QMC model 
is a natural extension of our calculation and work in this direction is under consideration. 
We also plan to investigate the effects of exchange terms for both the phase transition and 
the droplet properties. 

ACKNOWLEDGMENTS 

This work was supported in part by the University of Washington (USA), CNPq and 
FAPESP (Brazil) and FCT (Portugal) under the projects PRAX1S/2/2.1/F1S/451/94 and 
PRAXlS/P/FlS/12247/98. M.N. would like to thank the Centro de Fisica Teorica da Uni- 
versidade de Coimbra for its hospitality and financial support during her stay in Portu- 
gal. The authors also acknowledge useful discussions with Dr. C.A. Bertulani and Dr. 
M.B. Pinto. 



11 



REFERENCES 



[1] C.J. Pethick, D.G. Ravenhall, CP. Lorenz, Nucl. Phys. A 584 (1995) 675. 

[2] D.P. Menezes and C. Providencia, Nucl. Phys.A 650 (1999) 283. 

[3] D.P. Menezes and C. Providencia, Phys. Rev. C 60 (1999) 024313. 

[4] J. Boguta and A. R. Bodmer, Nucl. Phys. A 292 (1977) 413; A.R. Bodmer and C.E. 

Price, Nucl. Phys. A 505 (1989) 123; P.-G. Reinhard, M. Rufa, J. Maruhn, W. Greiner 

and J. Friedrich, Z. Phys. A 323 (1986) 13. 
[5] P.A.M. Guichon, Phys. Lett. B 200 (1988) 235. 
[6] B.D. Serot and J.D. Walecka, Adv. Nucl. Phys. 16 (1986) 1. 

[7] M. Arneodo, Phys. Rep. 240 (1994) 301; D.F. Geesaman, K. Saito and A.W. Thomas, 

Annu. Rev. Nucl. Part. Sci. 45 (1995) 337 . 
[8] I. Zakout and H. R. Jaqaman, Phys. Rev. C 59 (1999) 962. 

[9] S. Pal, M. Hanauske, I. Zakout, H. Stocker and W. Greiner, Phys. Rev. C 60 (1999) 
015802. 

[10] R.J. Furnstahl, B.D. Serot and H.-B. Tang, Nucl. Phys. A 598 (1996) 539; 615 (1997) 

441; 618 (1997) 446. 

[11] S. Fleck, W. Bentz, K. Shimizu and K. Yazaki, Nucl. Phys. A 510 (1990) 731. 
[12] K. Saito and A.W. Thomas, Phys. Lett. B 327 (1994) 9 . 

[13] P.A.M. Guichon, K. Saito, E. Rodionov and A.W. Thomas, Nucl. Phys. A 601 (1996) 
349. 

[14] P.G. Blunden and G.A. Miller, Phys. Rev. C 54 (1996) 359. 

[15] K. Saito, K. Tsushima and A.W. Thomas, Nucl. Phys. A 609 (1996) 339. 

[16] X. Jin and B.K. Jennings, Phys. Lett. B 374 (1996) 13; Phys. Rev. C 54 (1996) 1427; 

H. Miiller and B.K. Jennings, Nucl. Phys. A 640 (1988) 55. 

[17] P. K. Panda, A. Mishra, J. M. Eisenberg and W. Greiner, Phys.Rev. C 56 (1997) 3134; 

I. Zakout and H. R. Jaqaman, Phys.Rev. C 59 (1999) 968. 

[18] G. Krein, A.W. Thomas and K. Tsushima, Nucl. Phys. A 610 (1999) 313. 

[19] M.E. Bracco, G. Krein and M. Nielsen, Phys .Lett. B 432 (1998) 258. 

[20] A.W. Thomas, Nucl. Phys. A 629 (1998) 20; T. Frederico, B.V. Carlson, R.A. Rego, 

M.S. Hussein, J. Phys. G15 (1989) 397. 
[21] H. Reinhardt and H. Schulz, Nucl. Phys. A432 (1985) 630. 
[22] H. Miiller and R.M. Dreizler, Nucl. Phys. A 563 (1993) 649. 

[23] D. Von Eiff, J. M. Pearson, W. Stocker and M. K.Weigel, Phys. Lett. B 324 (1994) 279. 

[24] M. Centelles, M. Del Estal and X. Viiias, Nucl. Phys. A 635 (1998) 193. 

[25] P-G. Reinhard, Rep. Prog. Phys. 52 (1989) 439. 

[26] M. Barranco and J.R. Buchler, Phys. Rev. C 22 (1980) 1729. 

[27] H. Miiller and B.D. Serot, Phys. Rev. C 52 (1995) 2072. 

[28] B. D. Serot and J. D. Walecka, Phys. Lett. 87B (1979) 172. 

[29] M. Nielsen and J. Providencia, J. Phys. G 16 (1990) 649. 

[30] M. Centelles, X. Vifias, Nucl. Phys. A 563 (1993) 173. 

[31] M. M. Sharma, S. A. Moszkowski and P. Ring, Phys. Rev C 44 (1991) 2493. 
[32] G. Baym, H.A. Bethe and C.J. Pethick, Nucl. Phys. A 175 (1971) 225. 



12 



TABLES 



TABLE L Comparison of some properties of a droplet with total particle number ^4 = 40 
and proton fraction Np/A — 0.5 obtained in the Thomas-Fermi approximation with the 
QMC and NLl. The values of K, M* and are also given for both models. 



model 


Ppir = 0)fm"^ 


Ssurf (MeV/fm3) 


t (fm) 




K (MeV) 


M* (MeV) 


ms (MeV) 


QMC 


0.079 


1.15 


2.37 


4.00 


280 


754 


418 


NLl 


0.075 


1.08 


2.79 


3.94 


212 


535 


492.25 
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FIG. 1. Proton (lower curves) and neutron (upper curves) chemical potentials as a function of 
the proton fraction for P = 0.1 MeV/ fm^ for the QMC model (solid lines), for the Walecka model 
(dotted lines), and for the NLl (dashed lines). 
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FIG. 4. The droplet density surface energy as a function of the central proton fraction for the 

QMC model with (dashed line) and without (solid line) the electromagnetic field, and for the NLl 
without the electromagnetic field (dotted line). Results for a droplet with a total of 20 particles. 
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FIG. 5. The droplet density surface energy as a function of the central proton fraction for the 
QMC model with the electromagnetic field for a droplet with a total of 40 (solid line) and 20 
(dashed line) particles. 
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